Variation in the Antibacterial and Antioxidant Activities of Essential Oils of Five New Eucalyptus urophylla S.T. Blake Clones in Thailand.

Eucalyptus oils are widely used for a variety of purposes. This study investigates the terpenoid compositions and antibacterial and antioxidant activities of eucalypt leaf oils extracted from four E. urophylla clones and one E. urophylla × E. camaldulensis hybrid clone grown in Thailand. According to GC/MS analysis, the E. urophylla oils were mainly composed of 1,8-cineole, α-terpinyl acetate, ß-caryophyllene, and spathulenol, while 1,8-cineole, α-terpinyl acetate, p-cymene, and γ-terpinene were mostly identified in the hybrid oil. All eucalypt oils exhibited a significant bacteriostatic effect against Gram-positive bacteria, Streptococcus pyogenes, Staphylococcus aureus, Listeria monocytogenes, and Bacillus cereus. Only the hybrid oil had an effect on all Gram-negative bacteria tested, including Salmonella typhi, Escherichia coli, Pseudomonas aeruginosa, and Enterobacter aerogenes. These oils have antibacterial properties that vary according to their terpenoid content. Only the hybrid oil had a potent antioxidant effect, with an IC50 value of 4.21 ± 0.35 mg/mL for free radical (DPPH) scavenging. This oil's antioxidant effect may be a result of the phenolic terpenoids, thymol and carvacrol. As a result, these oils may be a novel source of antibacterial and antioxidant agents. Additionally, the antibacterial and antioxidant capabilities of the E. urophylla × E. camaldulensis hybrid essential oil are reported for the first time.

Changes in taxonomic and functional diversity of plants in a chronosequence of Eucalyptus grandis plantations.

Tree plantations have become one of the fastest-growing land uses and their impact on biodiversity was evaluated mainly at the taxonomic level. The aim of this study was to analyze environmental changes after the Eucalyptus plantation in an area originally covered by natural grasslands, taking into account the alpha and beta (taxonomic and functional) diversity of plant communities. We selected nine plantation ages, along a 12 years chronosequence, with three replicates per age and three protected grasslands as the original situation. At each replicate, we established three plots to measure plant species cover, diversity and environmental variables. Results showed that species richness, and all diversity indices, significantly declined with increasing plantation age. Canopy cover, soil pH, and leaf litter were the environmental drivers that drove the decrease in taxonomic and functional diversity of plants through the forest chronosequence. Based on the path analyses results, canopy cover had an indirect effect on plant functional diversity, mediated by leaf litter depth, soil pH, and plant species richness. The high dispersal potential, annual, barochorous, and zoochorous plant species were the functional traits more affected by the eucalypt plantations. We recommend two management practices: reducing forest densities to allow higher light input to the understory and, due to the fact that leaf litter was negatively associated with all diversity facets, we recommend reducing their accumulation or generate heterogeneity in its distribution to enhance biodiversity.

Population Divergence along a Genetic Line of Least Resistance in the Tree Species Eucalyptus globulus.

The evolutionary response to selection depends on the distribution of genetic variation in traits under selection within populations, as defined by the additive genetic variance-covariance matrix (G). The structure and evolutionary stability of G will thus influence the course of phenotypic evolution. However, there are few studies assessing the stability of G and its relationship with population divergence within foundation tree species. We compared the G-matrices of Mainland and Island population groups of the forest tree Eucalyptus globulus, and determined the extent to which population divergence aligned with within-population genetic (co)variation. Four key wood property traits exhibiting signals of divergent selection were studied-wood density, extractive content, and lignin content and composition. The comparison of G-matrices of the mainland and island populations indicated that the G-eigenstructure was relatively well preserved at an intra-specific level. Population divergence tended to occur along a major direction of genetic variation in G. The observed conservatism of G, the moderate evolutionary timescale, and close relationship between genetic architecture and population trajectories suggest that genetic constraints may have influenced the evolution and diversification of the E. globulus populations for the traits studied. However, alternative scenarios, including selection aligning genetic architecture and population divergence, are discussed.

Phylogenomics shows lignotuber state is taxonomically informative in closely related eucalypts.

Plant taxa can be broadly divided based on the mechanisms enabling persistence through whole-crown disturbances, specifically whether individuals resprout, populations reseed, or both or neither of these mechanisms are employed. At scales from species through to communities, the balance of disturbance-response types has major ramifications for ecological function and biodiversity conservation. In some lineages, morphologically identical populations except for differences in a disturbance-response trait (e.g. ±â€¯lignotuber) occur, offering the opportunity to apply genetic analyses to test whether trait state is representative of broader genetic distinctiveness, or alternatively, variation in response to local environmental conditions. In eucalypts, a globally-significant plant group, we apply dense taxon sampling and high-density, genome-wide markers to test monophyly and genetic divergence among pairs of essentially morphologically-identical taxa excepting lignotuber state. Taxa differing in lignotuber state formed discrete phylogenetic lineages. Obligate-seeders were monophyletic and strongly differentiated from each other and lignotuber-resprouters, but this was not the case for all lignotuber-resprouter taxa. One lignotuber state transition within our sample clade was supported, implying convergence of some non-lignotuber morphology characters. Greater evolutionary rate associated with the obligate-seeder disturbance-response strategy offers a plausible explanation for these genetic patterns. Lignotuber state is an important taxonomic character in eucalypts, with transitions in lignotuber state having contributed to the evolution of the exceptional diversity of eucalypts in south-western Australia. Differences in lignotuber state have evolved directionally with respect to environmental conditions.

GC-MS and HS-SPME-GC×GC-TOFMS Determination of the Volatile Composition of Essential Oils and Hydrosols (By-Products) from Four Eucalyptus Species Cultivated in Tuscany.

Essential oils are widely used as functional ingredients for potential multi-purpose functional uses. Hydrosols, co-products of the distillation of plant material, are used in food and cosmetic industries and in biological agriculture, but their volatile composition is poorly investigated. The volatile fractions of essential oils and hydrosols from four less-studied 1,8-cineol-rich Eucalyptus species (E. parvula L.A.S. Johnson & K.D. Hill, E. cinerea F. Muell, E. pulverulenta Sims and E. pulverulenta baby blue Sims), cultivated in Tuscany in a system of organic farming, were characterized by solvent dilution (essential oils) or extraction (hydrosols) followed by GC-MS and by HS-SPME-GC×GC-TOFMS analysis. GC-MS analysis showed that essential oils were mainly constituted by oxygenated monoterpenes, particularly 1,8-cineole, with monoterpenes hydrocarbons up to 10.8%. Relative differences in the abundance of minor terpenes as limonene, α-pinene, γ-terpinene, p-cymene, terpinen-4-ol, α-terpineol, and alloaromandrene were pointed out and seem to be suitable for differentiation among EOs of the four different Eucalyptus species. Hydrosols of these species were characterized for the first time: they were mainly constituted by oxygenated monoterpenes (97.6⁻98.9%), with 1,8-cineole up to 1.6 g/L, while monoterpene and sesquiterpene hydrocarbons were detected only in traces. HS-SPME-GC×GC-TOFMS analysis also allowed providing metabolic profiling of hydrosols for the direct comparison and visualization of volatile components, pointing out the potentially different uses of these products as functional ingredients in food, beverage, and cosmetic industries.

Impact of intense disturbance on the structure and composition of wet-eucalypt forests: A case study from the Tasmanian 2016 wildfires.

Fire is a key process in eucalypt communities, exerting a strong influence on the composition, structure and functioning of forests. Much of the research on the fire response of temperate, wet-sclerophyll trees in Australia comes from Victoria, where the dominant eucalypt is Eucalyptus regnans. In contrast, central and northern Tasmanian forests, dominated by Eucalyptus delegatensis, are relatively understudied. There is a need to determine whether Tasmanian wet-sclerophyll forests, though the same forest type in name, are functionally different in floristics and response to fire. Here we document the forest community response to a natural wildfire event in Tasmania-using opportunistic before/after control/impact (BACI) data from pre-existing monitoring plots. Uniting pre- and post-fire floristic data, we quantified mortality and regeneration of eucalypt, acacia and other dominant tree species, and tree ferns, Dicksonia antarctica, in response to wildfire. We also evaluated the density of eucalypt and acacia seedling establishment between burnt and unburnt forests, and quantified faunal responses to fire. Despite moderate-to-high intensity burning in patches across the plot, mortality of eucalypts, acacias and tree ferns due to fire were low. By contrast, fire-sensitive rainforest species showed low survival, though were able to persist in unburnt refugia. Eucalypt and acacia seedling regeneration was high in the burnt plot, suggesting that E. delegatensis forests regenerate without stand-replacing fire events. This contrasts to Victorian E. regnans forests, whose persistence is dependent on high-severity stand-replacing events. We also found some group-specific avifaunal and invertebrate responses to the fire event, which are broadly reflective of responses documented in other Victorian-based studies. Our results have implications for Tasmanian wet-forest silvicultural practices, which are based on the principle of stand-replacement after fire. The broader relevance of this work to forest ecology is in demonstrating the serendipitous opportunities that can arise with baseline monitoring plots.

An Array of Bioactive Compounds From Australian Eucalypts and Their Relevance in Pancreatic Cancer Therapeutics.

Pancreatic cancer (PC) is one of the most devastating human cancers, and despite the significant advances in the current therapeutic options, the overall survival rate for PC has remained static for the past 50 years. Plant-derived bioactive compounds play a vital role in cancer therapeutics by providing new lead compounds for future drug development. Therefore, the isolation, characterization, and identification of new bioactive compounds for the prevention and treatment of cancer continue to be an important aspect of natural product research. Many in vitro and in vivo studies published in the last few decades have established strong links between the phytochemical profile of eucalypts and anticancer activity. However, only a small number of these reports have attempted to demonstrate a relationship between the biological activity of eucalypt extracts and PC. This review focuses on potential anti-PC effects of an array of bioactive compounds present in various species of eucalypts. It also highlights the necessity for further in vitro and in vivo studies to develop a complete understanding of the potential this group of plants has for the development of potent and specific chemotherapeutic drugs for PC.

Speciation in the presence of gene flow: population genomics of closely related and diverging Eucalyptus species.

Speciation is a complex process that is fundamental to the origins of biological diversity. While there has been considerable progress in our understanding of speciation, there are still many unanswered questions, especially regarding barriers to gene flow in diverging populations. Eucalyptus is an appropriate system for investigating speciation mechanisms since it comprises species that are rapidly evolving across heterogeneous environments. We examined patterns of genetic variation within and among six closely related Eucalyptus species in subgenus Eucalyptus section Eucalyptus in south-eastern Australia (commonly known as the "green ashes"). We used reduced representation genome sequencing to genotype samples from populations across altitudinal and latitudinal gradients. We found one species, Eucalyptus cunninghamii, to be highly genetically differentiated from the others, and a population of mallees from Mount Banks to be genetically distinct and therefore likely to be a new undescribed species. Only modest levels of differentiation were found between all other species in the study. There was population structure within some species (e.g., E. obstans) corresponding to geographical factors, indicating that vicariance may have played a role in the evolution of the group. Overall, we found that lineages within the green ashes are differentiated to varying extents, from strongly diverged to much earlier stages of the speciation continuum. Furthermore, our results suggest the green ashes represent a group where a range of mechanisms (e.g., reticulate evolution and vicariance) have been operating in concert. These findings not only offer insights into recent speciation mechanisms in Eucalyptus, but also other species complexes.

Quantitative Trait Loci (QTLs) for Intumescence Severity in Eucalyptus globulus and Validation of QTL Detection Based on Phenotyping Using Open-Pollinated Families of a Mapping Population.

Intumescence is a nonpathogenic physiological disorder characterized by leaf blistering. This disorder can affect growth and development in glasshouses and growth chambers and may be confused with pathogenic diseases. We used quantitative trait loci (QTL) analysis to examine the genetic basis of variation in intumescence severity in Eucalyptus globulus, and test for colocation with previously detected QTLs for pathogen susceptibility. QTL analysis used the phenotype means of open-pollinated (OP) families of an outcrossed F2 mapping family (OP F3; n = 300) of E. globulus and the linkage map constructed in the F2. We validate this phenotyping approach for QTL analysis by assessing a trait previously used for QTL discovery in the F2 and showing the same major QTL was detected with the OP F3. For intumescence severity, five putative QTLs were detected across four linkage groups. Four of these did not colocate with previously reported QTLs for fungal pathogen susceptibility in Eucalyptus, suggesting the mechanisms underlying susceptibility to intumescence and to the two fungal pathogens are largely independent. This study demonstrates there is a genetic basis for variation in intumescence severity, reports the first QTL for intumescence severity in plants, and provides a robust framework for investigating the potential mechanisms involved.

Seedling response to environmental variability: The relationship between phenotypic plasticity and evolutionary history in closely related Eucalyptus species.

PREMISE OF THE STUDY: Phenotypic plasticity is an important means through which organisms cope with environmental variability. We investigated seedling plasticity in the green ash eucalypts within a phylogenetic framework to examine the relationship between plasticity and evolutionary history. The green ashes are a diverse group, which include the tallest flowering plant in the world (Eucalyptus regnans) and a rare mallee less than 1 m tall (E. cunninghamii). METHODS: Seedlings of 12 species were exposed to high and low nutrient and water availability in a factorial experiment. Leaf trait and total plant plasticity were evaluated using the phenotypic plasticity index. A phylogeny of the species was estimated using genome-wide scans. KEY RESULTS: We found significant differences in functional traits across species, growth forms, and substrates in response to changes in resource availability. Many traits (e.g., leaf width) were highly plastic for most species. Interspecific differences in leaf-level plasticity was significant, however plasticity was not correlated with phylogeny. Species with broader environmental niches had higher leaf-level plasticity than species with narrower environmental ranges. CONCLUSIONS: Plastic responses to environmental variability can differ widely among closely related species, and plasticity is therefore likely to be associated with many factors, including habitat and range size, as well as evolutionary history. Our results provided insights for species delimitation in Eucalyptus, which have management implications. Because of the high number of rare species and that other species are commercially important, a more comprehensive understanding of plasticity is essential for predicting their response to changing climates.

Genetic divergence among Psidium accessions based on single nucleotide polymorphisms developed for Eucalyptus.

The goal of this study was to analyze the genetic divergence among Psidium species accessions based on SNPs developed for Eucalyptus. Fifty-three Psidium accessions, including 47 P. guajava, were genotyped with EUCHIP60K. The dendrogram similarity ranged from 0.58 to 1.00, with a cophenetic value of 0.97. Five groups were identified at dendrogram cut point of 0.7: the first with 44 guava accessions, the second with 1 guava accession, the third with 3 P. guineense accessions, the forth with 2 guava accessions, and the fifth with 3 P. cattleianum accessions. The Bayesian analyses suggested seven subpopulations, with formation of two additional groups with guava accessions. Primers designed with Eucalyptus SNP sequences resulted in reliable Psidium amplicons on 6% polyacrylamide gels. In general, the SNP dendrogram agreed with biological genus structure, since different species were not grouped, indicating that transferability among Myrtaceae genus was possible and reliable.

Genomic Scans across Three Eucalypts Suggest that Adaptation to Aridity is a Genome-Wide Phenomenon.

Widespread species spanning strong environmental (e.g., climatic) gradients frequently display morphological and physiological adaptations to local conditions. Some adaptations are common to different species that occupy similar environments. However, the genomic architecture underlying such convergent traits may not be the same between species. Using genomic data from previous studies of three widespread eucalypt species that grow along rainfall gradients in southern Australia, our probabilistic approach provides evidence that adaptation to aridity is a genome-wide phenomenon, likely to involve multiple and diverse genes, gene families and regulatory regions that affect a multitude of complex genetic and biochemical processes.

Climate and geochemistry as drivers of eucalypt diversification in Australia.

Eucalypts cover most of Australia. Here, we investigate the relative contribution of climate and geochemistry to the distribution and diversity of eucalypts. Using geostatistics, we estimate major element concentrations, pH, and electrical conductivity at sites where eucalypts have been recorded. We compare the median predicted geochemistry and reported substrate for individual species that appear associated with extreme conditions; this provides a partial evaluation of the predictions. We generate a site-by-species matrix by aggregating observations to the centroids of 100-km-wide grid cells, calculate diversity indices, and use numerical ecology methods (ordination, variation partitioning) to investigate the ecology of eucalypts and their response to climatic and geochemical gradients. We find that ß-diversity coincides with variations in climatic and geochemical patterns. Climate and geochemistry together account for less than half of the variation in eucalypt species assemblages across Australia but for greater than 80% in areas of high species richness. Climate is more important than geochemistry in explaining eucalypts species distribution and change in assemblages across Australia as a whole but there are correlations between the two sets of environmental variables. Many individual eucalypt species and entire taxonomic sections (Aromatica, Longistylus of subgenus Eucalyptus, Dumaria, and Liberivalvae of subgenus Symphyomyrtus) have distributions affected strongly by geochemistry. We conclude that eucalypt diversity is driven by steep geochemical gradients that have arisen as climate patterns have fluctuated over Australia over the Cenozoic, generally aridifying since the Miocene. The diversification of eucalypts across Australia is thus an excellent example of co-evolution of landscapes and biota in space and time and challenges accepted notions of macroecology.

Insights into temperature modulation of the Eucalyptus globulus and Eucalyptus grandis antioxidant and lignification subproteomes.

Eucalyptus grandis and Eucalyptus globulus are among the most widely cultivated trees, differing in lignin composition and plantation areas, as E. grandis is mostly cultivated in tropical regions while E. globulus is preferred in temperate areas. As temperature is a key modulator in plant metabolism, a large-scale proteome analysis was carried out to investigate changes in the antioxidant system and the lignification metabolism in plantlets grown at different temperatures. Our strategy allowed the identification of 3111 stem proteins. A total of 103 antioxidant proteins were detected in the stems of both species. Hierarchical clustering revealed that alterations in the antioxidant proteins are more prominent when Eucalyptus seedlings were exposed to high temperature and that the superoxide isoforms coded by the gene Eucgr.B03930 are the most abundant antioxidant enzymes induced by thermal stimulus. Regarding the lignin biosynthesis, our proteomics approach resulted in the identification of 13 of the 17 core proteins involved in this metabolism, corroborating with gene predictions and the proposed lignin toolbox. Quantitative analyses revealed significant differences in 8 protein isoforms, including the ferulate 5-hydroxylase isoform F5H1, a key enzyme in catalyzing the synthesis of sinapyl alcohol, and the cinnamyl alcohol dehydrogenase isoform CAD2, the last enzyme in monolignol biosynthesis. Data are available via ProteomeXchange with identifier PXD005743.

Analysis of cellulose synthase gene expression strategies in higher plants using RNA-sequencing data.

The transcriptomes from different organs and tissues of western poplar, eucalyptus, soybean and common bean were studied. The expression level of cellulose synthase genes was notably different in different types of tissues and organs in studied plants. For common bean and eucalyptus transcriptome the domination of certain cellulose synthase genes was typical. These prevailing genes made up more than 50 % of the total expression pull of cellulose synthases. On the contrary, cellulose synthase expression pulls of wes-tern poplar and soybean were distributed between multiple genes. The different expression strategies of CesA-genes may reflect a phylogenetic processes that occurred in genomes studied.

Vessel diameter and related hydraulic traits of 31 Eucalyptus species arrayed along a gradient of water availability.

We present two comprehensive data sets that describe xylem vessel diameters and related sapwood traits for species of Eucalyptus from arid and semi-arid woodlands and forests in Australia. Between 2009 and 2014, sapwood of mature trees was sampled in south-western, south-eastern and eastern Australia. One additional species was sampled from tropical north-western Australia. The first data set describes samples collected from the basal stem section (130 cm above ground) of three individuals of 31 species of which eight species were replicated at sites that differed in climatic conditions. The second data set describes vessel characteristics of three trees from each of 10 species that were sampled at 8 m below the tree apex. The sampled trees of these 10 species are also part of the first data set. In total, we report diameters (D) for over 25 100 vessels, from 494 digital images taken from 117 trees. We also report vessel frequencies, void-to-wood ratios, sapwood densities and hydraulically weighted vessel diameters (Dh). Supporting data of the first data set include tree diameter at breast height (130 cm above ground), tree height, sample locations, and summary climate data. In this data set, diameter of individual vessels ranges from 10 to over 300 µm, and vessel frequency from 360 to 9070 vessels cm-2 . Wood density ranges from 0.47 to 0.96 g cm-3 . Void-to-wood ratio ranges from 6% to 27% and Dh ranges from 46 to 236 µm. Mean annual rainfall (P) at sample sites ranges from 246 to 2274 mm and FAO56 reference evaporation (E) from 777 to 2110 mm. The aridity index (E/P) ranges from 0.15 to 2.93 (dimensionless). Tree diameters range from 9 to 90 cm and tree heights range from 6 to 70 m. D and Dh in the second data set range from 11 to 271 and 68 to 205 µm, respectively. These datasets will make a valuable contribution to future continental-scale and global-scale studies of the relationship between xylem hydraulic architecture and climate. The data sets are unique in the sense that they are phylogenetically constrained, allowing in-depth assessment of plasticity of hydraulic attributes within a single tree genus.

High density, genome-wide markers and intra-specific replication yield an unprecedented phylogenetic reconstruction of a globally significant, speciose lineage of Eucalyptus.

We used genome-wide markers and an unprecedented scale of sampling to construct a phylogeny for a globally significant Eucalyptus lineage that has been impacted by hybridisation, recent radiation and morphological convergence. Our approach, using 3109 DArT markers distributed throughout the genome and 540 samples covering 185 terminal taxa in sections Maidenaria, Exsertaria, Latoangulatae and related smaller sections, with multiple geographically widespread samples per terminal taxon, produced a phylogeny that largely matched the morphological treatment of sections, though sections Exsertaria and Latoangulatae were polyphyletic. At lower levels there were numerous inconsistencies between the morphological treatment and the molecular phylogeny, and taxa within the three main sections were generally not monophyletic at the series (at least 62% polyphyly) or species (at least 52% polyphyly) level. Some of the discrepancies appear to be the result of morphological convergence or misclassifications, and we propose some taxonomic reassessments to address this. However, many inconsistencies appear to be the products of incomplete speciation and/or hybridisation. Our analysis represents a significant advance on previous phylogenies of these important eucalypt sections (which have mainly used single samples to represent each species), thus providing a robust phylogenetic framework for evolutionary and ecological studies.

A simple protocol to determine lignin S/G ratio in plants by UHPLC-MS.

A simple extraction protocol and an ultra-high performance liquid chromatography-mass spectrometry (UHPLC-MS) method for the determination of the syringyl/guaiacyl (S/G) ratio in lignin is reported herein. The method was entirely developed using stems of three Eucalyptus species, which were hydrolyzed with NaOH and partitioned with ethyl ether; vanillin (from the G monomer) and syringaldehyde (from S monomer) were quantified. The S/G ratios obtained were comparable to those usually reported for eucalyptus. The data for one of the eucalyptus species were compared with those obtained with a widely accepted method using thioacidolysis and gas chromatography-mass spectrometry (GC-MS). The method was also applied to sugarcane and showed to be reliable. The yield of the NaOH hydrolysis of the monolignols ranged from 89.94 to 95.69%, with more than 77.12% of recuperation in the liquid-liquid extraction. The whole analytical procedure was validated, achieving results with less than 4.38% of variation. The lowest LOD and LOQ were 0.01 and 0.05 µg/mL, respectively. In addition, the method combines reliability and a fast and direct quantification.

The Eucalyptus terpene synthase gene family.

BACKGROUND: Terpenoids are abundant in the foliage of Eucalyptus, providing the characteristic smell as well as being valuable economically and influencing ecological interactions. Quantitative and qualitative inter- and intra- specific variation of terpenes is common in eucalypts. RESULTS: The genome sequences of Eucalyptus grandis and E. globulus were mined for terpene synthase genes (TPS) and compared to other plant species. We investigated the relative expression of TPS in seven plant tissues and functionally characterized five TPS genes from E. grandis. Compared to other sequenced plant genomes, Eucalyptus grandis has the largest number of putative functional TPS genes of any sequenced plant. We discovered 113 and 106 putative functional TPS genes in E. grandis and E. globulus, respectively. All but one TPS from E. grandis were expressed in at least one of seven plant tissues examined. Genomic clusters of up to 20 genes were identified. Many TPS are expressed in tissues other than leaves which invites a re-evaluation of the function of terpenes in Eucalyptus. CONCLUSIONS: Our data indicate that terpenes in Eucalyptus may play a wider role in biotic and abiotic interactions than previously thought. Tissue specific expression is common and the possibility of stress induction needs further investigation. Phylogenetic comparison of the two investigated Eucalyptus species gives insight about recent evolution of different clades within the TPS gene family. While the majority of TPS genes occur in orthologous pairs some clades show evidence of recent gene duplication, as well as loss of function.

Patterns of Reproductive Isolation in Eucalyptus-A Phylogenetic Perspective.

We assess phylogenetic patterns of hybridization in the speciose, ecologically and economically important genus Eucalyptus, in order to better understand the evolution of reproductive isolation. Eucalyptus globulus pollen was applied to 99 eucalypt species, mainly from the large commercially important subgenus, Symphyomyrtus. In the 64 species that produce seeds, hybrid compatibility was assessed at two stages, hybrid-production (at approximately 1 month) and hybrid-survival (at 9 months), and compared with phylogenies based on 8,350 genome-wide DArT (diversity arrays technology) markers. Model fitting was used to assess the relationship between compatibility and genetic distance, and whether or not the strength of incompatibility "snowballs" with divergence. There was a decline in compatibility with increasing genetic distance between species. Hybridization was common within two closely related clades (one including E. globulus), but rare between E. globulus and species in two phylogenetically distant clades. Of three alternative models tested (linear, slowdown, and snowball), we found consistent support for a snowball model, indicating that the strength of incompatibility accelerates relative to genetic distance. Although we can only speculate about the genetic basis of this pattern, it is consistent with a Dobzhansky-Muller-model prediction that incompatibilities should snowball with divergence due to negative epistasis. Different rates of compatibility decline in the hybrid-production and hybrid-survival measures suggest that early-acting postmating barriers developed first and are stronger than later-acting barriers. We estimated that complete reproductive isolation can take up to 21-31 My in Eucalyptus. Practical implications for hybrid eucalypt breeding and genetic risk assessment in Australia are discussed.